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Abstract 

We present analytical solutions to the nonlinear equations describing 
the behavior of a gas of neutrinos with two flavors. Self-maintained 
coherent flavor oscillations are shown to occur when the gas density 
exceeds a critical value determined by the neutrino masses and the mean 
neutrino energy in the gas. Similar oscillations may have occurred in 
the early Universe. 



To appear in Physical Review D, July 1995 



*Permanent address: Physics Department, City College of New York, New York, NY 10031, USA. 
E-mail : samuel@scisun.sci . ccny. cuny. edu 



I. Introduction 

The properties of neutrinos passing through matter have attracted much attention 
in the past few years. The reason is that the large electron-lepton number of matter 
results in flavor-dependent neutrino propagation and, at suitable densities, enhanced 
neutrino flavor oscillation The effect can be attributed primarily to neutrino- 
electron forward scattering through W ± exchange. Among the consequences is a 
possible resolution of the solar neutrino problem 0. 

Another situation with lepton imbalance is the early Universe at a tempera- 
ture below 1 MeV, where neutrino propagation is affected by the excess of electrons 
over positrons. Some implications for the early Universe of neutrino oscillations en- 
hanced by W ± exchange have been considered in M. In this scenario, however, 
the self-interactions of neutrinos through Z° exchange cannot be neglected because 
the neutrino density is relatively high. Indeed, numerical simulations of the early 
Universe show that neutrino properties can be significantly modified by these effects 
P, |6|, 0| . In the parameter region of 0, there is even a novel neutrino-flavor conversion 
mechanism that is different from the MSW effect. 

The neutrino self- interactions are nonlinear. In a general context, they must be 
studied numerically. Moreover, in the early Universe, neutrino behavior is controlled 
by a combination of factors. In addition to the electron-positron imbalance and the 
presence of the neutrino gas itself, these include other effects such as the expansion 
rate of the Universe. In this paper, we eliminate these additional complications by 
considering a simplified situation consisting of a homogeneous gas of self-interacting 
neutrinos in a box of fixed volume V, with no other leptons present. Our goal is to 
obtain results that are analytical but that nonetheless describe nonlinear features of 
neutrino behavior. 

We consider here two situations: the pure neutrino gas (Sec. 3), and a gas con- 
taining both neutrinos and antineutrinos (Sec. 4). Throughout the paper, we assume 
that hard-scattering processes are negligible compared to forward scattering. This is 
valid provided the energy E of the neutrino satisfies the condition GfE 2 /(he) 3 <C 1. 
Forward scattering corresponds to phase-interference effects and hence to neutrino 
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oscillations. Under these conditions, neutrinos are neither created nor destroyed but 
are simply transformed from one flavor to another. 

For simplicity in what follows, we restrict ourselves to oscillations between electron 
and muon neutrinos. For vacuum oscillations, relevant parameters are the vacuum 
mixing angle 9 and the mass-squared difference A = (m| — m 2 ) c 4 . The effective 
hamiltonian describing free neutrino propagation is diagonal in the mass-eigenstate 
basis. However, weak-interaction processes produce and destroy flavor-eigenstate 
neutrinos. The mass-eigenstate neutrinos v\ and u 2 are related to the left-handed 
flavor-eigenstate neutrinos v e L and u^l by 

v \ — u eL cos 9 — v^l sin 9 , v 2 = v eL sin 9 + cos 9 , (1.1) 

where mi and rri2 are the masses of v\ and v 2 , respectively. 

The vacuum-oscillation period T A is given by 

47rEh , . 

T A = ^- . (1.2) 

For a gas of N v neutrinos with a finite energy spread, neutrinos oscillate with var- 
ious periods. In the limit of negligible interactions, a neutrino (or mixed neutrino- 
antineutrino) gas with large N u exhibits oscillation decoherence, i.e., the net neutrino- 
flavor content is time independent at late times. This constant asymptotic behavior 
is independent of initial conditions. For example, suppose one begins with a gas 
of electron neutrinos. Some will convert into muon neutrinos so that time-varying 
behavior occurs initially. The flavor content of an individual neutrino at time t is 
1 -sin 2 29 (1 - cos (2nt/T A )) for u e and sin 2 29 (1 - cos (2nt/T A )) for The summa- 
tion over cosine functions with various periods Ta leads asymptotically to a constant 
function of t if sufficiently many terms are present, i.e., if N v is sufficiently large. 
The ratio of electron neutrinos to muon neutrinos in the limit t — > oo becomes the 
constant factor (1 — sin 2 29)/ sin 2 29. 

Vacuum behavior dominates provided the neutrino gas is sufficiently dilute, so that 
interactions due to Z Q exchange remain unimportant. The dominance is controlled 
by a dimensionless parameter k, defined by 

A 



K 



2V2G F En v 



2 



where n u = N u /V is the neutrino density. This parameter can be written as the ratio 
T u /T Al where 

T v = -= (1.3) 

is the time scale associated with neutrino interactions. When the neutrino density is 
low, k is large. Neutrino-neutrino forward scattering occurs infrequently compared 
to a vacuum oscillation period, and the behavior is similar to a non-interacting gas. 
This region of parameter space is characterized by decoherence. 

In contrast, when the neutrino density is large so that k <C 1, neutrino interac- 
tions are important. Many neutrino-neutrino interactions occur during a vacuum- 
oscillation period. Numerical simulations for the pure neutrino gas || reveal the exis- 
tence in this parameter region of a collective mode of the nonlinear dynamics in which 
the behaviors of individual neutrinos are correlated. Significant numbers of neutri- 
nos oscillate in unison. We refer to this counterintuitive behavior as self-maintained 
coherence. Self-maintained coherence is also seen in numerical simulations of a gas of 
neutrinos and antineutrinos |5|, ||, |7j . A system consisting initially of electron neutri- 
nos does not decohere. Instead, oscillatory behavior is observed, even at late times. 
A primary goal of this paper is to obtain an analytical description of self-maintained 
coherence for k <§C 1. 

If m2 < m i so that A < (or, alternatively, if A > and 9 > 7r/2), then there is a 
large region in the A-8 parameter space for which self-maintained coherence emerges 
for neutrino oscillations in the early Universe ||. The behavior begins smoothly, but, 
due to the expansion of the Universe and the varying electron and positron densities, 
coherent oscillations emerge at around 100 seconds after the big bang. Self-maintained 
coherent oscillations may thus have played a role in early-Universe physics. 

Throughout the rest of this paper, we work in units with h = c = 1 . 

II. Background Material 

A single relativistic neutrino oscillating in vacuum obeys the equation 

dv 

i- = Hv , (2.1) 
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where u(t) is the two- component flavor wave function 

"-(£) < 2 - 2 ' 

with z/*z/ e + = 1, and where the effective hamiltonian H is given by 

ml + ml /l 0\ A /-cos 20 sin20\ 

4£ VO 1 / 4E V sin20 cos 20/ " 1 dJ 

The probability for the particle to be an electron neutrino is v*v e , while that to be a 
muon neutrino is v*v... 

A convenient and standard vector reformulation of the above equations exists 
It is useful both for visualization and for numerical simulation of oscillations. Define 
the vector 

v = (W - ufr, 2Re(X*^), 2Im(^)) . (2.4) 



Then, the neutrino oscillation equation ( |2.3| ) is equivalent to that governing a particle 
of unit mass and charge moving in a magnetic field B given by 

S = tE • < 2 ' 8 ' 

where 

A = A (cos 29, - sin 29, 0) . (2.6) 



For an antineutrino, Eqs. ( |2 . 3| ) — ( P76|) hold if neutrino wave functions are replaced by 
antineutrino wave functions, i.e., v — ► v. We denote the corresponding vector for an 
antineutrino by w. Throughout this work, we use the expression "magnetic field" to 
refer to an effective magnetic field, as opposed to a physical one. 

For a gas, there is a vector v? for the j'th neutrino and a vector w k for the fcth 
antineutrino. The equations governing the self-interacting gas become || [Kj] 

HiP 

— = iP x Bl (2.7) 
for vectors associated with neutrinos, and 

^ = ™ k xBl (2.8) 
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for antineutrinos. Here, the magnetic fields B 3 V and B^ are given by 

Bi = ^-V„, St = ^ + VX , (2.9) 

k 

where the energy of the fcth antineutrino is denoted E . An asterisk on a vector 
indicates a change in sign of the third component. The vacuum contributions to the 
magnetic fields are the terms in Eq. ( |2.9| ) dependent on A, given in Eq. ( |2.6| ). The 
potential V vv is generated by Z° exchange and is given by 

= ((v) - (w*)) , (2.10) 

where Gf — 1-17 x 1CT 11 MeV -2 is the Fermi coupling constant, and 

(rJ> = Xy, (w) = J2^ ■ (2-11) 

3 k 

— * _ _ 

In the absence of V w , the first-order differential equations in (|2.7|) and ( j2.8|) de- 
couple and are linear. The system is then solvable and the solution corresponds to 
a non-interacting gas in which each neutrino undergoes vacuum oscillatory behavior. 
When V uu is present, the equations in ( |2.7|) and Q2.8|) are both coupled and nonlinear. 
For this reason, we refer to V vv as the nonlinear term. This nonlinearity leads to 
interesting effects. 

Individual neutrinos and antineutrinos are neither created nor destroyed under 
our assumptions. The equations expressing this, 

d (v j ■ v j ) d (w k -w k ) 



follow from Eqs. (|2.7|) and 

Different normalizations of iP are possible. Above, we have chosen v 3 ■ v 3 = 1 
along with the interpretation that the index j labels individual neutrinos. For this 
case, the index j ranges from 1 to N v , where N v is the total number of neutrinos. 
A second normalization convention follows from noting that neutrinos with the same 
energy obey the same oscillation equation. One can therefore perform a sum over 
all vectors of the same energy. With this second normalization, \v 3 \ represents the 
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number of neutrinos of energy EK The index j then ranges over the possible energy 
values. Similar normalization choices exist for antineutrinos. 

For both the above normalization conventions, the total number of neutrinos N v 
and antineutrinos N 9 is given by 

N U = Y.\V\, N P = J2\w k \ • (2.13) 

j k 

The neutrino and antineutrino densities n v and n p can then be obtained by dividing?] 
by V, that is, n u = N u /V, n D = N p /V. 

For numerical purposes the second normalization scheme is more useful. For the 
mathematical treatment in the current work, we use the first convention with tP -tP = 1 
and w k ■ w k = 1 . 

III. The Pure Neutrino Gas 

In this section, we analyze self-maintained coherence for a pure neutrino gas. For 
this system, N D = and the contribution to the neutrino potential becomes 

?~ = ^<8> • (3-D 

For defmiteness, we consider the situation in which N u electron neutrinos are placed 
in the box at time t = 0, so that the initial conditions are 

^'(0) = (1,0,0) . (3.2) 

At t = 0, the ratio k- 7 of the vacuum term to the neutrino-neutrino term for the jth 
neutrino is 

= _ . (3.3) 

2y/2G F Ein v 

When k? <C 1, the vacuum term is dominated by the neutrino- neutrino interaction 
term, and self-maintained coherence appears in computer simulations 0. Neglecting 
the non-linear term, neutrinos with larger energies oscillate slower and neutrinos with 



1 Both normalization schemes discussed here can also be modified by multiplying neutrino vectors 
by 1/V. In this situation, neutrino vectors become densities, and V in the nonlinear term V vv given 
in Eq. ( |2.1C| ) must be replaced by 1. 
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smaller energies oscillate faster. However, a large V vv term boosts slow neutrinos and 
retards fast neutrinos. 

To obtain an analytical solution, we can take advantage of a feature of the motion 
called alignment JF]: numerical simulation shows that vectors in the nonlinear system 
point in a common direction when the k? are small. Alignment of the jth neutrino 
implies that the approximation 

]f ( 3 - 4 ) 

is good. Here, f v (t) is the average neutrino vector. This feature suggests we seek an 
analytical solution for r v (t). 

An equation for r v (t) can be obtained by summing over j in Eq. ( p.7|) and using 
Eqs. (|2J), and (gT 



dr v A 

— — = r\, x , 3.5) 

dt 2E ' K ' 

where the average inverse energy 1/E is defined by 

— = — V-L . (3.6) 

Equation ( |3.5|) shows that the self-maintained coherence in the pure neutrino gas 
is formally equivalent to the oscillation of a single neutrino in vacuum. Hence, the 
average neutrino vector undergoes vacuum oscillations with an effective energy Eq. 

It is useful for later purposes, when the neutrino-antineutrino gas is considered, 
to display the solution of Eq. ( |3.5|) . The initial conditions for r v are 

f%(0) = (1,0,0) . (3.7) 

Equation ( |3.5| ), when written in components, is 

dt \2E J d ' 

dr v2 ( A \ 

cos 20 r v3 , 



dt \2Eo 

= -— - (r,,i sin 26* + r„ 2 cos 26») . (3. 
dt 2E 
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These equations simplify in the vacuum-mass-eigenstate basis denoted by R(t) and 
given by 

Ri = n cos 29 - r 2 sin 29 , R 2 = n sin 29 + r 2 cos 29 , R 3 = r 3 . (3.9) 



In this basis, the equations resemble those in Eq. (|3.8|) with 8 = 0. Thus, Ri(t) is a 
constant. The equations for R 2 and R% combine to give a harmonic-oscillator system. 
Incorporating the initial conditions (|3.7|), we find 



Ri(t) = cos 26 , 
A 



R 2 (t) = sin 26 cos 



2En 



\2E 

Returning to the flavor basis, we obtain the desired solution: 



R 3 (t) = - sin 20 sin (—t) . (3.10) 



ri(t) = cos 2 26 + sin 2 26 cos ( ——t 

2E 

A 



r 2{t) = — sin 26 cos 26 ^1 — cos 



2K 



o 



r 3 (t) = -sin 20 sin ( — - 1 ) . (3.11) 



2E Q 

Summarizing, the solution in the dense neutrino parameter region is given by Eqs. 
!4|) and (|3.11|) . These equations describe self-maintained oscillations. All neutrinos 



oscillate in unison. Note that perfect alignment is obtained in the limit kP — ► 0. 

When some & are large, the corresponding neutrinos do not participate in the 
collective mode. If most neutrinos have small & then self-maintained coherence still 
occurs but with a smaller amplitude. The criterion for self-maintained coherence for 
a pure neutrino gas is Ko < 1, where 

K = . 3.12 

2V2G F E n u v ; 

We have compared our analytical solution to numerical simulations. Excellent 
agreement is obtained when all h? 1. Even for the case in Figure 8 of ref. ||, for 
which 10% of the neutrinos had k j > 1, agreement between the analytical approach 
and numerical simulations is to about 5% for the oscillation period and the amplitude. 
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An intuitive understanding of alignment and self-maintained coherence is as fol- 
lows. Assume that most neutrino vectors are aligned. These vectors point along the 
average vector f(t) and collectively rotate around A. Consider a particular neutrino 
with a higher energy than average. Let v be its vector. In the absence of the non- 
linear term, v rotates around A at a relatively slow rate. Suppose v begins to lag f. 
Then, because the nonlinear term is much bigger than the vacuum term, the neutrino 
experiences a large magnetic field in the direction of f(t). Consequently, v rotates 
around r. After half a period, v will have rotated to a position leading the group. 
Hence, v cannot lag behind or otherwise separate from the group. 

A similar argument holds for any neutrino with energy lower than average. If 
v begins to lead r, it experiences a large magnetic field in the direction of f(t) and 
so rotates around r rather than A. It follows that neutrinos with energies different 
from the average do not rotate around A at varying rates but instead stay together 
in a group. This is alignment. Since the group follows r, which rotates around A, 
oscillatory behavior arises. This is self-maintained coherence. 

IV. The Neutrino- Antineutrino System 

In this section, we study self-maintained coherence for a dense gas containing 
both neutrinos and antineutrinos. Numerical simulations reveal that alignment holds 
separately for neutrinos and antineutrinos || Consequently, the approximations 

= ^(t)-^^(t) (4.1) 

for the j'th neutrino and the fcth antineutrino are good. By summing over j and k in 
Eqs. (|2.7]) and (|2.8| ), differential equations for r v and r w are obtained: 

f^ x (i-^l ■ < 4 - 2 > 

f^„x(J, + C] , (4.3) 



where 



V w = V2G F (n v r v - n„r*) . (4.4) 
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In Eq. ( |4.2| ), 1/E is the average inverse neutrino energy ( |3.6| ), and 1/E is the 
analogous quantity for antineutrinos. 

For definiteness, we consider the situation with an equal number of electron neu- 
trinos and antineutrinos placed in the box at time t = 0, n v = n p . Then, the initial 
conditions are 

^(0) = (1,0,0) , w k (0) = (1,0,0) . (4.5) 

Also, for simplicity we take the antineutrinos to have the same average inverse energy 
as neutrinos: Eq = Eq. This holds, for example, in the more restricted case when 
the energy distributions of neutrinos and antineutrinos are the same, E = E J for all 
k = j. 

The symmetry of the initial conditions suggests the ansatz 

r w {t)=r v {t) . (4.6) 

It follows that Eqs. (|4.2|) and ( |4.3| ) are equivalent, so it suffices to solve one of the 
pair to demonstrate consistency of the ansatz. The nonlinear term in Eq. ( |4 . 4| j only 
has a third component, V uv = 2\f2n v Gpr v3 (0,0, 1). The problem therefore reduces 
to solving the equations 

dr v \ ( A 



dt \2E 
dr v2 ( A 



dt \2E 
dr v3 A 



sin 28 j r v3 - 2\f2n v G F r v2 r v3 , 

cos 26j r v3 + 2V2n v G F r 

r v i sin 26 -\- r v2 cos 2$) (4.7) 



dt 2E 

determining the components of average neutrino vector. 

These equations again simplify in the vacuum-mass-eigenstate basis R(t) given in 
Eq. ( |3.9| ). The vector R obeys Eq. ( |4.7| ) for r with 8^0. In what follows, it is 
convenient to make the further change of variables 

R^t) = yi (s) , R 2 (t) = y 2 (s) , R 3 (t) = ^y 3 (s) , (4.8) 

where k is given in Eq. (|3.12|) and where 



" ( • " = lvi • (4 - 9) 
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The oscillation equations simplify to 



dy 



,- = -y2V3 , (4.io) 

ds 

dy2 kq 

-ds~ = ym + Y m ' (4 - n) 

— = -y 2 . (4.12) 
ds 

The initial conditions ( |4.5| ) become 

R(Q) = y(0) = (cos 29, sin 29, 0) . (4.13) 

To proceed, we can take advantage of the conservation of neutrino number, which 
in the present variables is expressed as 

yl + yl + ^vi = i • (4.14) 

This equation, which is a consequence of Eqs. CN= . 1 0|) — (N= • 1 3|) , specifies y 2 in terms of 
yi and y^. Furthermore, an equation for y% in terms of y\ is obtained by substituting 
Q4.12j ) into CN= - 1 0Q and integrating. The above observations determine y 2 and y% in 
terms of y\ as 



y 2 = ±Vl-yi-«o(yi- cos 29) , 2/3 = ±V 2 (m-cos20) . (4.15) 

Specifying the signs corresponds to specifying different stages of the motion, as dis- 
cussed below (see Eq. ( |4.23| ) ) . 

At this point, we need only obtain y\(t). A differential equation for this variable 
can be found by differentiating ( |4.10| ) with respect to s, using Eqs. ( |4.11| ) and ( |4.12| ), 



incorporating Eqs. ( f4.14j ) and ( |4 . 1 5| ) , multiplying by dyi/ds, and integrating once. 



The integration constant is fixed using Eq. ( 4. 13| ) . This procedure gives 



1 / dyi 

2 \ ds 



2 

I =-V(yi) , (4.16) 



where 

- V (i/i) = (yi - cos 29) (l - y\ - k ( Vl - cos 29)) . (4.17) 
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Equation ( 4.16Q is analogous to one describing the motion of a particle of unit mass 



moving in the potential V (yi) associated to an anharmonic oscillator with a cubic 
term. 

The motion of y\ is between y™ m and y™ ax , where 

yf n = cos 26 , j/™* = x + + cos 26 . (4.18) 

The point x + determines one of the two zeros of the potential — V(x + cos20) = 
— x (x — x + ) (x + xJ), where 



x ± = ^1 + /cocos20 + k§/4=fcos20=fk o /2 • (4.19) 

In the region y™ n < y\ < y™ ax , it follows that — V (yi) > 0, so both sides of Eq. 
( |4.16| ) are positive. 



The solution of Eq. (|4.16|) is by quadrature in x. Changing the integration variable 
using x = x + w 2 gives the implicit solution 




x_ Jo 



2 f\ 7~ aw 



(4.20) 



- w 2 ) (1 + q 2 w 2 ) 

where q 2 = x + /x_. This expression can be inverted for y±(s) using the sine-amplitude 
and delta-amplitude Jacobi elliptic functions sn and dn, defined as 

f sn(u,k) fa fdn(v,k) fa 

U = I , , V = 







1 - 1 2 ) (1 - k 2 t 2 ) J o ^J(l-t 2 )(l-(l-k 2 )t 2 ) 

(4.21) 

We finally obtain 



sin 2 



(26) sn 2 (Ms, 

yi (s)=cos26 + — [ZU) 9 ) V A V A\ . ( 4 - 22 ) 



where xq = x + + x_ = 2y 1 + cos 20 + Kq/4. 
The motion consists of four stages: 

stage 1 ■ ^ n ^ n -> n g& <r n 

stage 2 



stage 3 
stage 4 



y 2 >o, y 3 <0, ^>0, ^<0, ^<0, 

Z/ 2 <0, y 3 <0, ^<0, ^<0, §>0, 

Z/ 2 <0, y 3 >0, ^>0, ^>0, ^>0, 

2/ 2 >0, y 3 >0, ^<0, %>0, ^<0 



(4.23) 
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Since y\ > cos 29, it follows that y\ > for all stages of the motion. During the 
motion, y™ n < y 2 < y™ ax , where y™ ax = sin 29 = —y™ m . 

The motion of y is roughly circular about A. Recall that A points along the 1-axis 
of the vacuum-mass-eigenstate basis. At t — 0, y— (cos 29, sin 29, 0) = {y 1 \ m , y™ '*, 0), 
so y points in the direction of the 1-axis of the flavor basis. During stage 1, y drops 
below the 1-2 plane. It then passes to stage 2 when it drops below A. At the 1- 
to-2 transition point, yi = y™ ax and yi = 0. The vector y continues its motion 
below the 1-2 plane during stage 2, until y\ and 2/2 obtain their minimum values and 
ys = 0. In stages 3 and 4, the motion is reversed, except that y% is positive so that 
y is above the 1-2 plane. The maximum value of yi is again achieved at the 3-to-4 
transition, where y is above A and y<i = 0. The cycle is completed when y returns to 
y = (cos 29, sin 29, 0). Hence, one entire cycle of motion involves two cycles of y\. The 
signs of 1/2 and 2/3 in Eqs. ([4.15| ) are determined from the second and third columns 
in Eq. flO^). 

Figures la-c display the behavior over two periods of each of the three components 
of the vector r for the case with sin 2 29 = 0.81 and = 0.1. The time scale is plotted 
in terms of the s variable, which is equivalent to the time t measured in units of l//i. 
It is convenient to use s so that A and E need not be specified (compare with Eq. 
( [4.9[ )). Extra oscillations appear in the second component, plotted in Fig. lb. They 
arise from the projection of the orbit onto the 2 axis. The variables yi, y<i, y% have 
no such effects. Figure 2 shows the same orbit in a three-dimensional plot. 

In general, the half-period of y± in s is determined from Eq. (|4.20|) by setting 
'2/1 — y\ ax ■ Since a complete cycle involves two y\ cycles, one obtains 

- dw (4.24) 



for the period S vD in s. Equation (fO|) then implies that the period T v9 in t is 



T up = ^ . (4.25) 
Hence, the period T vD is of the order of the geometric mean of the time scales as- 



sociated with the vacuum and nonlinear terms: T vV ~ ^JT^T V . The motion for the 
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neutrino-antineutrino gas is thus on the order of 1/ <J~Kq times faster than in the pure 
neutrino case of Sect. 3. 

Anther interesting feature of the behavior of the neutrino-antineutrino system is 



its near planarity, apparent since R% is related to 2/3 by a factor of ^Jk /2. The range 
of -R3 is of order ^/k^ <C 1. Hence, the bulk of the motion is in the 1-2 plane. Planarity 
is a feature observed in the numerical simulations of refs. |§, §], • In the variables R 
or r, the orbit is similar to a highly eccentric ellipse. 

V. Summary 

In this paper, we have provided analytical solutions to the nonlinear equations 
describing the behavior of a gases containing two flavors of neutrinos, both with and 
without antineutrinos. For a dense pure neutrino gas, the solution is given by Eqs. 
( |3.4|) and (|3.11 ), while for a dense neutrino-antineutrino gas the solution is specified 



by Eqs. (p|), (pD, fj), (14151) , O and (pp . 



The behavior of the neutrino-antineutrino gas differs from that of the pure neutrino 
case. The former is controlled by elliptic functions, whereas the latter is governed by 
trigonometric functions. Our analytical results agree in detail with prior numerical 
simulations in the region with k < 1. 

We have demonstrated analytically that self-maintained coherent flavor oscilla- 
tions occur when the gas density exceeds a critical value, given in terms of the mean 
neutrino energy and the neutrino masses. Oscillations of this type may have occurred 
in the early Universe. 
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FIGURE CAPTIONS 

Figure 1. Components of the vector r as a function of scaled time, s = fit, for the 
case sin 2 29 = 0.81 and k = 0.1. (a) The component r\. (b) The component r 2 . (c) 
The component r 3 . 

Figure 2. The three-dimensional orbit for the case sin 2 26 = 0.81 and k = 0.1. 
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